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Abstract. Entanglement represents a pure quantum effect involving two or more particles. 
Spin systems are good candidates for studying this effect and its relation with other collective 
phenomena ruled by quantum mechanics. While the presence of entangled states can be easily 
verified, the quantitative estimate of this property is still under investigation. One of the most 
useful tool in this framework is the concurrence whose definition, albeit limited to S = 1/2 
systems, can be related to the correlators. We consider quantum spin systems defined along 
chains and square lattices, and described by Heisenberg-like Hamiltonians: our goal is to clarify 
the relation between entanglement and quantum phase transitions, as well as that between the 
concurrence the and the specific quantum state of the system. 



1. Introduction 

The occurrence of collective behavior in many-body quantum systems is asso- 
ciated with classical and quantum correlations. The latter, whose name is entan- 
glement, cannot be accounted for in terms of classical physics and represents the 
impossibility of giving a local description of a many-body quantum state. Entan- 
glement is expected to play an essential role at quantum phase transitions (QPT), 
where quantum effects manifest themselves at all length scales, and the problem 
has recently attracted an increasing interest [1, 2, 3, 4, 5, 6]. 
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Moreover, entanglement comes into play in quantum computation and commu- 
nication theory, being the main physical resource needed for their specific tasks [7]. 
In this respect, the perspective of manipulating entanglement by tunable quantum 
many-body effects appears intriguing. 

In this paper, entanglement estimators are found to give important insight in 
the physics of S = 1/2 spin systems, for which the concurrence give quantitative 
definition of bipartite entanglement. Quantum spin chains and two dimensional 
lattices in external fields are studied. Two striking features are found: the occur- 
rence of a factorized ground state at a field hf and that of a QPT at h c > hf , 
where multipartite entanglement plays an essential role. 

The entanglement estimators have been calculated by numerical simulations, 
carried out in models with linear dimension L , via Stochastic Series Expansion 
(SSE) Quantum Monte Carlo based on a modified directed-loop algorithm [8]. We 
have verified that the inverse temperature j3 = 2L is suitable to test the T = 
behaviour. 



We consider the antiferromagnetic S =1/2 XYZ model in a uniform magnetic 



where J>0 is the exchange coupling, (ij) runs over the pairs of nearest neigh- 
bors, and h=g^BH/J is the reduced magnetic field. The canonical transformation 
Si' y ^(-iySi' y with I = 1(2) for i belonging to sublattice 1(2), transforms the 
coupling in the xy plane from antiferromagnetic to ferromagnetic. The Hamilto- 
nian (1) is the most general one for an anisotropic S = 1/2 system with exchange 
spin-spin interactions. However, as real compounds usually display axial symme- 
try, we will henceforth consider either A J/ =1 or A z =l. Moreover, we will apply 
the field along the z-axis, i.e. h=(0, 0, h). 

This paper focuses on the less investigated case A a =l, defining the XYX model 
in a field. Due to the non-commutativity of the Zeeman and the exchange term, for 
A^t^I this model is expected to show a field-induced QPT on any D-dimensional 
bipartite lattice, with the universality class of the D-dimensional Ising model in 
a transverse field [9]. The two cases A y <l and Ay>l correspond to an easy- 
plane (EP) and easy-axis (EA) behavior, respectively. The ordered phase in the 
EP(EA) case arises by spontaneous symmetry breaking along the x(y) direction, 
which corresponds to a finite value of the order parameter M x (M y ) below the 
critical field h c . At the transition, long-range correlations are destroyed, and 
the system is left in a partially polarized state with field-induced magnetization 
reaching saturation only as h—too. This picture has been verified so far in D=l 



2. 



The model 



field: 
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only, both analytically [10]and numerically [2, 11]. 

Besides its quantum critical behavior, a striking feature of the model Eq. (1) is 
the occurrence of an exactly factorized ground state for a field hf(A y ) lower than 
the critical field h c . This feature was previously predicted in Ref. [12] for magnetic 
chains, and its entanglement behaviour has been studied in Ref. [2]. Our QMC 
simulations have given evidence [3] for a factorized ground state to occur also in 
D=2. We have then rigorously generalized the proof of factorization to the most 
general Hamiltonian Eq. (1) on any 2D bipartite lattice. The proof will be soon 
reported elsewhere, but we here outline the essential findings in the following. The 
occurrence of a factorized ground state is particularly surprising if one considers 
that we are dealing with the S = 1/2 case, characterized by the most pronounced 
effects of quantum fluctuations. However, in the class of models here considered, 
they are fully uncorrelated[12] at h = /if, thus leading to a classical-like ground 
state. 



3. Entanglement estimators 

In order to calculate the entanglement of formation [13] in the quantum spin 
system described by Eq. (1) we make use of the one-tangle and of the concurrence. 
The one-tangle [14, 15] quantifies the T = entanglement of a single spin with the 
rest of the system. It is defined as: 

n = 4detpW, = (I + Y,M a a a )/2; (2) 

a 

is the one-site reduced density matrix, M a = (S a ), a a are the Pauli matrices, 
and a = x,y, z. In terms of the spin expectation values M a , one has: 

r 1 = l-4^(M Q ) 2 . (3) 

a 

The concurrence [16] quantifies instead the pairwise entanglement between two 
spins at sites i and j, both at zero and finite temperature. For the model of 
interest, in absence of spontaneous symmetry breaking (M x = 0) the concurrence 
reads [15] 

dj = 2max{0,C;.,C;;.} , (4) 

where 

Cfii = »U \ - «u <1?J • (5) 

c£- = IC + 4"l-/(i + 4*) 2 -( M2 ) 2 > 

and gfj* = (S^S") are the static correlators. 
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One-tangle and concurrence are related by the Coffman-Kundu-Wootters (CKW) 
conjecture [14], recently proved by Osborne and Verstraete [17], stating that 

Tl>T 2 =£Cg, (7) 

which expresses the crucial fact that pairwise entanglement does not exhaust the 
global entanglement of the system, as entanglement can also be stored in 3-spin 
correlations, 4-spin correlations, and so on. Therefore n-spin entanglement and 
m-spin entanglement with m^n are mutually exclusive and this is really a unique 
feature of entanglement as a form of correlation. The difference with classical 
correlations is evident. 

Due to the CKW conjecture, the entanglement ratio [2] R=T2/t\<\ quantifies 
the relative weight of pairwise entanglement, and its deviation from unity shows the 
relevance of n-spin entanglement with n>2. Although indirect, the entanglement 
ratio is a powerful tool to estimate multi-spin entanglement. 



4. The factorized state 

In the one dimensional system described by Eq. (1) the occurence of an exactly 
factorized ground state for a field hf(A y ) lower than the critical field h c , was 
predicted in Ref. [12]. In the case of the XYX model, this factorizing field is 

hf = yj2(l + A y ) . (8) 

At h = hf the ground state of the model takes a product form |^) = <S)fLi \i>i)i 
where the single-spin states are eigenstates of (n^) ■ S), where n 1 ( 2 ) being the 
local spin orientation on sublattice 1 (2). Taking n = (cos <f> sin 8, sin (f> sin 9, cos 6), 
one obtains [12] 



01 = O,0 2 = ^, 0i = 6 2 = cos" 1 yj(l + A y )/2 . (9) 

In particular, for h = h{ the spin orientation in the quantum state is exactly 
the same as in the classical limit of the model made of continuous spins with 
effective spin length S e a = 1/2; this means that quantum fluctuations only set the 
length of the effective classical spin. The factorized state of the anisotropic model 
continuously connects with the fully polarized state of the isotropic model in a 
field for A y = 1 and h = 2. 

In the two-dimensional case, we have found that for any value of the anisotropies 

and A z , there exists an ellypsoid in field space 

hi ^ ''" -+ 7 a !f\ ^T = 4 (10) 



(l + Aj / )(l + A Z ) (l + Ay)(Ay + A Z ) (l + A z )(A y + A z ) 
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such that, when h lies on its surface, the ground state of the corresponding model 
is factorized, \*&)= ®iLi |Y>i). The single-spin states are eigenstates of (nj-S), 
n/ being the local spin orientation on sublattice /. We will hereafter indicate with 
hf (factorizing field) the field satisfying Eq. (10); at h=hf, the reduced energy 
per site is found to be e=— (l+A^+A^/^. In the particular case of A z =l and 

h=(0, 0,/i), the factorizing field takes the simple expression hf =2^2(1+ A^). As 
for the structure of the ground state, taking n/=(cos 4>i sin 9j, sin sin 9j, cos 9i), 
the analytical expressions for 0/ and Oi are available via the solution of a system 
of linear equations. 

The local spin orientation turns out to be different in the EP and EA cases, 
being 

01=0 , 02=VT , #1 = 02= COS" 1 y/(l + Ay)/2 , (11) 

for A y <l, and 

0i=vr/2 , 02=-vr/2 , 1= 6 2 = cos" 1 ^2/(l+A y ) , (12) 

for A y >l. 

Despite its exceptional character, the occurrence of a factorized state is not 
marked by any particular anomaly in the experimentally measurable thermody- 
namic quantities. On the other hand, entanglement estimators are able to pin 
down the occurrence of such factorized states with high accuracy, as shown in the 
following section. 

5. Entanglement and quantum phase transitions 

The entanglement estimators display a marked anomaly at the factorizing field, 
where they clearly vanish, as expected. When the field is increased above hf, the 
ground-state entanglement has a very steep recovery, accompanied by the QPT at 
h c > hf. The system realizes therefore an interesting entanglement effect controlled 
by the magnetic field. This feature is associated with the many-body behavior of 
the system and it is shown in Fig.l and Fig.2 for one and two-dimensional systems, 
respectively. 

The concurrence terms, are generally short-ranged, and usually do not extend 
farther than the third neighbor. The longest range of Cij is indeed observed around 
the factorizing field hf [18, 19]. 

The sum of squared concurrences T2 is seen to be always smaller than or equal to 
the one-tangle n, both for one- and two-dimensional systems. This is in agreement 
with the CKW conjecture [14]. The total entanglement is only partially stored in 
two-spin correlations, and it appears also at the level of three-spin entanglement, 
four-spin entanglement, etc. In particular, we interpret the entanglement ratio 



R = t 2 /ti , 



(13) 




Fig. 1: One dimensional XYX model with A y = 0.25 : t\ and T2 versus h. Insets: 
entanglement ratio T2/T1 close to the QPT. . 



as a measure of the fraction of the total entanglement stored in pairwise cor- 
relations. This ratio is plotted as a function of the field in the insets of both 
figures. A striking anomaly occurs at the quantum critical field h c , where R dis- 
plays a very narrow dip. According to our interpretation, this result shows that 
the weight of pairwise entanglement decreases dramatically at the quantum critical 
point in favour of multi-spin entanglement. In fact, due to the CKW conjecture, 
and unlike classical-like correlations, entanglement shows the special property of 
monogamy, namely full entanglement between two partners implies the absence 
of entanglement with the rest of the system. Therefore multi-spin entanglement 
appears as the only possible quantum counterpart to long-range spin-spin corre- 
lations occurring at a QPT. This also explains the somewhat puzzling result that 
the concurrence remains short-ranged at a QPT while the spin-spin correlators 
become long-ranged [1], and it evidences the serious limitations of concurrence as 
an estimate of entanglement at a quantum critical point. In turn, we propose the 
minimum of the entanglement ratio R as a novel estimator of QPT, fully based on 
entanglement quantifiers. 




Fig. 2: Two dimensional XYX model: t\ and t% versus h. Left panel: easy-plane 
case (A y = 0.25); Right panel: easy-axis case (A y = 4). Insets: entanglement 
ratio tz/tx close to the QPT.. 

6. From spin configurations to entanglement properties 

We here analyze the entanglement of formation between two spins, as quantified 
by the concurrence C, in terms of spin configurations. In the simplest case of two 
isolated spins in the pure state \4>) the concurrence may be written as C = \ J2i a ?l> 
where a« are the coefficients entering the decomposition of \<j)) upon the magic 
basis. In this case, it can be easily shown [18] that C extracts the information 
about the entanglement between two selected spins by combining probabilities 
and phases relative to some specific two-spin states. 

In general, one should notice that a finite probability for two spins to be in 
a maximally entangled (Bell) state does not guarantee per se the existence of 
entanglement between them, since this probability may be finite even if the two 
spins are in a separable state. In a system with decaying correlations, at infinite 
separation all probabilities associated to Bell states attain the value of 1/4, but 
this of course tells nothing about the entanglement between them, which is clearly 
vanishing. It is therefore expected that differences between such probabilities, 
rather than the probabilities themselves give insight in the presence or absence of 
entanglement. 

When the many-body case is tackled, the mixed-state concurrence of the se- 
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lected spin pair has an involved definition in terms of the reduced two-spin density 
matrix [16]. We here assume that TC is real, has parity symmetry (meaning that 
either 7i leaves the z component of the total magnetic moment unchanged, or 
changes it in steps of 2), and is further characterized by translational and site- 
inversion invar iance. 

Let us select two specific spins in the system. We indicate by pg, with (3 = 
1,2,3,4, the probabilities for the two spins to be in one of the Bell states (maxi- 
mally entangled), where [3 = 1, 2 (3,4) refers to the parallel (antiparallel) ones. We 
do also introduce p l , p E , p m , p w , which represent the probabilities for the two spins 
to be in the (factorized) states of the standard base, where indexes I, II (III, IV) 
refers to the parallel (antiparallel) ones. 

It can be shown [18] that C' and C" (as defined in Eqs. (5) and (6) with 
indexes dropped for simplicity) can be written in terms of the above probabilities: 

2C ' = \p3-p4\-2VPA > (14) 
2C" = \pi - p 2 \ ~ (1 ~ Pi ~ P2) = 

= \Pl ~P2\- Z^PwPtv > ( 15 ) 

where we have used p m =p w and hence P3 + P4 = 2p w = 2^/p w p w . The expression 
for C" may be written in the particularly simple form 

2C" = 2max{pi,p2}-l , (16) 

telling us that, in order for C" to be positive, it must be either p\ > 1/2 or 
P2 > 1/2. This means that one of the two parallel Bell states needs to saturate 
at least half of the probability, which implies that it is by far the state where the 
spin pair is most likely to be found. 

Despite the apparently similar structure of Eqs. (14) and (15), understanding 
C' is more involved, due to the fact that ^Jp Y p^ cannot be further simplified unless 
Pi = Pn- The marked difference between C and C" reflects the different mecha- 
nism through which parallel and antiparallel entanglement is generated when time 
reversal symmetry is broken; (p : ^ p n and hence M z ^ 0). In fact, in the zero 
magnetization case, it is p n = p l = (p\+p2)/1 and hence 

2C" = 2max{p 3 ,p 4 } - 1 , (17) 

which is fully analogous to Eq. (16), so that the above analysis can be repeated 
by simply replacing p\ and p 2 with p^, and p±. 

For M z / 0, the structure of Eq. (17) is somehow kept by introducing the 
quantity 

so that 

2C' = 2max{p 3 ,P4} - (1 - A 2 ) , (19) 
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meaning that the presence of a magnetic field favors bipartite entanglement associ- 
ated to antiparallel Bell states, \e^) and In fact, when time reversal symmetry 
is broken the concurrence can be finite even if p$, p4 < 1/2. 

^From Eqs. (16) and (19) one can conclude that, depending on C being finite 
due to C or C", the entanglement of formation originates from finite probabilities 
for the two selected spins to be parallel or antiparallel, respectively. In this sense 
we speak about parallel and antiparallel entanglement. 

Moreover, from Eqs. (14) and (15) we notice that, in order for parallel (an- 
tiparallel) entanglement to be present in the system, the probabilities for the two 
parallel (antiparallel) Bell states must be not only finite but also different from 
each other. Thus, the maximally entangled states result mutually exclusive in the 
formation of entanglement between two spins in the system, which is in fact present 
only if one of them is more probable than the others. The case p\ = P2 = 1/2 
(P3 = Pa = 1/2) corresponds in turn to an incoherent mixture of the two parallel 
(antiparallel) Bell states. 

The above analysis suggests the first term in C" (C) to distill, out of all possible 
parallel (antiparallel) spin configurations, those which are specifically related with 
entangled parallel (antiparallel) states. These characteristics reinforce the meaning 
of parallel and antiparallel entanglement. 

7. Conclusions 

We have analyzed the behaviour of one- and two-dimensional 5 = 1/2 antifer- 
romagnets displaying a field-driven quantum phase transition. We have shown that 
while bipartite entanglement does not show any peculiar feature at the transition, 
the entanglement ratio R [2], which measure the relevance of bipartite entangle- 
ment with respect to the total entanglement content of the system, has a marked 
dip at criticality: this indicates that multipartite entanglement rules the QPT, R 
being a powerful tool to detect it. 

On the other hand, bipartite entanglement is found to efficiently detect classical- 
like ground states, even the highly non-trivial ones which are invisible under an 
analysis based upon standard magnetic observables. 

Finally, an interpretation of the analytical expression of the concurrence is 
given in terms of spin configurations, leading to a deeper insight into the relation 
between entanglement properties and state configurations in many-body systems. 
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